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SUBSUB  -  fl  Submarine  Engagement  Model 

James  N.  Eagle 

Department  of  Operations  Research 
Naval  Postgraduate  School 

Monterey,  CR  939*43 

Introduction 

This  document  reports  the  results  of  an  examination  of  the  submarine 
engagement  model  SUBSUB,  (uritten  by  Presearch,  Inc  of  Fairfax,  Vfi,  and 
modified  for  use  at  the  Naval  lUar  College  (NUJC),  Newport,  RI .  This 
model  considers  a  submarine  versus  submarine  interaction  from  search 
through  detection,  classification,  localization,  closure,  attack,  and 
counterattack.  The  submarine  missions  modelled  in  are  barrier  patrol, 
area  patrol,  and  transit.  The  version  of  SUBSUB  examined  runs  on  an 
IBM-PC.  Other  versions  exist  which  operate  on  larger  computers. 

Purpose  of  this  Study 

The  primary  purpose  of  the  examination  was  to  assess  the  modelling 
methods  used  in  SUBSUB  and  to  judge  the  program's  usefulness  for 
seminar  war  gaming  at  NUJC.  This  evaluation  was  conducted  using 
documentation  and  FORTRRN  source  code  supplied  by  Presearch,  Inc. 

R  secondary  purpose  of  the  examination  was  to  describe  the  basic 
models  involved  and  suggest  possible  improvements  or  extensions. 

Limits  of  this  Study 

SUBSUB  is  an  extensive  model,  or  more  accurately,  an  integrated 
collection  of  models.  The  examination  considered  in  detail  the  available 
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documentation  and  FORTRAN  source  code  for  the  search,  detection,  and 
closure  modules.  The  environmental  and  platform  data  bases  mere  not 
examined. 


SUBSU8  Oueruiem 

SUBSUB  models  one-on-one  barrier,  area,  and  transit  search  missions 
for  submarine  platforms.  It  is  an  analytical  model,  as  opposed  to  a 
discrete  time  step  simulation.  The  submarine  can  be  assigned  a  variety 
of  current  or  projected  sonar  systems  and  rneapon  capabilities.  Each 
platform  is  modelled  in  either  an  aggressive  or  evasive  posture 
Countermeasures  are  allomed.  Each  encounter  begins  uuith  detection 
and  continues  through  classification,  localization,  closure,  attack,  and 
counterattack.  Figure-of-merit  inputs  (target  source  levels  and 
environmental  parameters)  are  dramn  from  a  disk-based  data  base.  The 
version  of  SUBSUB  examined  contained  only  minter  environments,  but 
others  could  presumably  be  added  as  necessary. 

The  program  begins  mith  the  user  revieiuing  the  last  selections  for 
environment,  scenario,  platforms,  sensors,  and  meapons.  Any  or  all  of 
these  selections  can  be  changed.  The  possible  scenarios  for  the 
opposing  submarines  (called  Blue  and  Red)  are:  barrier  search/transit, 
area  search/transit,  area  search/area  search,  transit/transit,  and  direct 
support/penetrator. 

When  the  user  is  satisfied  mith  the  initial  conditions,  the  program 
calculations  are  begun.  Outputs  of  the  program  include  all  input  data, 
nominal  detection  and  counterdetection  ranges  (apparently  the 
maximum  range  at  mhich  mean  signal  excess  is  0),  and  a  sequence  of 
calculated  probabilities  for  Blue  and  Red  detection,  closure,  attack, 
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counterattack,  and  kill.  Random  variables  uniformly  distributed 
between  0  and  1  are  compared  to  each  calculated  probability.  If  the 
random  variable  is  less  than  the  probability,  the  event  in  question  is  said 
to  occur.  The  model  then  progresses  to  the  next  event  in  the 
engagement  event  tree.  For  example,  if  it  is  determined  that  Blue 
achieves  a  secure  detection,  the  next  probability  calculated  is  that  of 
correct  classification  by  Blue,  followed  by  closure,  localization,  and 
attack. 


Detailed  Model  Description 

Probability  Density  for  Detection  Range.  For  the  selected  location  and 
season,  a  propagation  loss  curve  is  retrieved  from  a  disk-based 
database.  Using  the  passive  sonar  equation  a  mean  figure-of-merit 
(FOM)  is  calculated.  The  actual  figure-of-merit  is  assumed  to  be 

FOM  =  FOM  +  i, 

where  £  is  a  normally  distributed  random  variable  with  mean  0  and 
standard  deviation  specified  by  the  program.  Using  the  selected 
propagation  Loss  curve,  a  detection  range  is  associated  with  each  value 
of  FOM.  In  this  manner  a  density  function  for  detection  range  is 
determined.  (Figure  1.) 


The  same  calculation  is  done  for  the  target,  resulting  in  a  density 
function  for  counterdetection  range.  Although  it  is  not  clear  from  the 
documentation  or  source  code,  the  range  associated  uuith  each  FOM  is 
most  likely  the  maximum  range  such  that  propagation  loss  equals  FOM. 

An  implicit  assumption  of  this  model  is  that  detection  and 
counterdetection  ranges  are  fixed  during  the  entire  encounter  at  values 
determined  by  probability  distributions.  This  is  in  contrast  to  discrete 
time  step  computer  simulations  using  lambda-sigma  or  Gauss-Markov 
detection  models,  tuhere  time-varying  fluctuations  in  the  environment 
cause  these  ranges  to  change  over  time. 

Distribution  of  Minimum  Range  Achieved  during  Search.  If  p(t)  is  the 
minimum  range  achieved  during  the  search  from  time  0  to  time  t,  then 
PCPfl(t,R)  is  defined  as 
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Prob{p(t)  x<  R}. 

The  functional  form  of  PCPfl  changes  with  the  geometry  of  the  encounter. 
In  the  barrier  search  problem,  PCPfl  is  assumed  only  a  function  of  R.  That 
is,  PCPfl(R)  is  the  probability  that  the  minimum  range  achieved  during  the 
barrier  penetration  is  R  or  less.  The  calculation  of  Pcpfl  for  barrier  and 
area  search  geometries  will  be  discussed  in  detail  later. 

Probability  of  Secure  Detection  (PCD).  The  probability  of  secure 

detection  is  computed  by  conditioning  on  the  detection  range  of  the 
searcher.  Specifically,  if  fRo(r)  is  the  density  function  for  detection  range 

and  Rcd  is  counterdetection  range,  then  the  probability  of  achieving  a 
secure  detection  by  time  t  is 

00 

Prob{RCDx<r}  fRD(r)  PCpp(t,r)  dr.  (1) 

o 

In  SUBSUB,  DCD  is  evaluated  by  conditioning  on  the  searcher  FOM. 

Letting  r(o>)  be  the  detection  range  associated  with  FOM  v  and  fFOM(o>)  be 
the  probability  density  function  for  FOM,  then  DCD  becomes 

CO 

Prob{RCD^r(^)}  fF0MM  PCPfl(t,rM)  du  (2) 

<J 

-00 

This  expression  is  integrated  numerically.  The  probability  of  a  secure 
counterdetection  is  computed  similarly. 

This  method  of  calculating  the  probability  of  a  secure  detection  has 
much  to  recommend  it.  Rs  opposed  to  the  procedures  of  reference  [  1  ], 
this  method  incorporates  the  relative  motion  between  the  searcher  and 
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target.  Thus  the  probability  of  achieving  a  secure  detection  in  an  area 
search  mission  mill  in  general  be  different  from  that  in  a  barrier  mission. 
This  occurs  because  the  distribution  of  CPR  ranges  is  different  for  the 
two  cases. 

The  disadvantage  of  this  approach,  however,  is  the  requirement  to 
perform  a  numerical  integration.  This  slows  the  execution  of  the 
program  But  it  should  be  noted  that  even  the  simpler,  geometry- 
independent  method  of  reference  (1]  requires  numerical  integration  (or 
extensive  table  look  up).  It  appears  that  this  is  the  numerical  price  that 
must  be  paid  for  secure  detection  calculations  when  detection  and 
counterdetection  ranges  are  random  variables. 

If  searcher  has  a  significant  acoustic  advantage,  then  these 
calculations  are  probably  not  necessary.  In  this  case,  the  probability  of 
secure  detection  by  time  t  approximately  equals  the  probability  of 
detection  by  time  t.  But  when  detection  and  counterdetection 
capabilities  are  nearly  equal,  then  it  becomes  more  important  to 
carefully  consider  the  target's  potential  to  counterdetect. 

Barrier  search  PrPR.  From  the  SUBSUB  documentation,  PCPfi(R)  for  a 
searcher  conducting  a  barrier  patrol  against  a  target  with  a  constant 
course  and  speed  is 


[2R/IBUU  sin  |3))  -  (R2/[BLLI  T  tan2|3))  +  [2R,  /(BUJ  T  tan  |3)) 


13] 


inhere 

fl1  =  ((90-p)/360)  7TR2  if  R  <  (BUJ-T]/2 
or 

fl1  =  5((BUU-T)/2]2  tan  cr  +  (90-|3-cr)7TR2/360  if  (BLU-T)/2  <  R  <  (BLU-T)/(2sin  (3) 
or 

fl1  =  .5[(BW-T)/2)2  (1/tan  (3)  if  R  >  (BLU-T)/[2sin  |3). 

The  other  definitions  are: 

T  =  luidth  of  searcher's  track, 

BUJ  =  barrier  luidth  (the  target's  penetration  is  uniformly  distributed 
over  this  length], 

Vq  =  searcher  speed, 

VT  =  target  speed, 
p  =  tan  1  (Vj/Vq), 
cr  =  cos'1  ((BULhT)/(2R)). 

This  expression  is  a  modification  of  the  "standard"  calculations  in 
Koopman's  reference  [2],  The  general  computational  method  used  here 
is  to  consider  the  barrier  penetration  in  "target-stationary  relative 
space".  That  is,  the  target  is  assumed  to  be  stationary,  and  all  relative 
speed  for  the  encounter  is  provided  by  the  searcher  The  searcher's 
relative  speed  component  across  the  barrier  front  is  the  searcher's 
actual  speed,V0.  Rnd  the  component  of  relative  speed  perpendicular  to 
the  barrier  is  the  target  speed,  VT.  Then  the  probability  that  the  target 
comes  Luithin  range  R  of  the  searcher  during  a  barrier  penetration  is  the 
ratio  of  the  area  "covered"  by  the  searcher  in  relative  space  during  one 
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pass  across  the  barrier  to  the  total  area  that  could  be  occupied  by  the 
target.  In  Figure  2.,  this  covered  area  is  shaded. 
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BW 

- ► 


Figure  2.  Calculation  of  PCPA  for  Barrier  Search. 
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There  are  several  problems  uuith  (3).  The  first  is  that  (3)  approaches 
negative  infinity  as  VT  becomes  small  or  R  becomes  large.  The  physical 
explanation  for  this  is  that  the  method  used  to  calculate  the  area  of  the 
shaded  region  of  Figure  2.  (i.e.,  dividing  the  region  into  parallelograms, 
triangles  and  sections  of  circles)  fails  for  small  VT  or  large  R.  The 
mathematical  reason  for  this  behavior  is  that  the  second  additive  term 
in  (3)  dominates  the  other  tuuo  terms  in  these  limiting  cases. 


By  way  of  example.  Figure  3.  is  a  plot  of  PCPfl,  as  calculated  by  (3),  for 
target  speeds,  VT,  of  1  to  10  knots.  Here  R  is  18  nautical  miles  (nm),  V0  is 
10  knots,  BUJ  is  100  nm,  and  T  is  60  nm.  Since  the  actual  PCPfi  should  be  a 
decreasing  function  of  VT,  the  calculated  values  for  VT  less  than  about 
2.5  knots  are  suspect. 
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Figure  4.  is  a  plot  of  PCPfl  vs.  range,  R.  In  this  example,  VT  is  5  knots;  R 
varies  from  20  to  60  nm;  and  V0,  BLU,  and  T  are  as  before,  flctual  PCPfl 
should  be  an  increasing  function  of  R  (as  long  as  PCPfl  <  1,  at  least)  but  is 
seen  here  to  decrease  for  R  greater  than  about  36  nm. 

Although  not  mentioned  in  the  supplied  documentation,  the  SUBSUB 
code  (in  subroutine  RCPR)  modifies  (3)  so  as  to  reduce  the  problem 
created  by  large  values  of  R.  Specifically,  luhen  R>(BUU-T)/(2  sin  p),  then 
Pcpfiis  calculated  tuith  R  set  to  min{R,(Ttan2p)/sinp}.  This  change 
prevents  PCPfl  from  approaching  negative  infinity,  but  can  introduce  a 
discontinuity  into  the  function  and  does  not  prevent  PCPfl  from 
decreasing  uuith  increasing  R. 
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Another  problem  with  the  approach  used  to  derive  (3)  is  that  it  is  not 
applicable  for  BUI  <  T.  In  fact,  the  BUJ  terms  in  the  denominators  of  the 
additive  terms  of  (3)  cause  the  calculated  Pcpfl  to  become  infinite  as  BLU 
approaches  0.  The  SUBSUB  program  does  not  warn  the  user  that  BLU  <  T 
is  not  allowed. 

A  third  (and  relatively  minor)  problem  with  (3)  is  that  if  the 
methodology  use  to  compute  A1  in  the  first  two  cases  is  followed  for  the 
third  case  (i.e.,  when  A  > (BUJ-T)/(2sin  p)),  then  A1  should  be 

.5{  ((BUI-D/2)2  (1/tan  p)  -  (tan  p)(R  -  (BUI-T)/(2sin  (3))2  }. 

However,  as  discussed  above,  using  this  methodology  at  all  is  not 
recommended. 

Alternative  analytical  barrier  models  are  available  which  avoid 
these  problems  and  are  more  general  in  their  application.  A  good 
example  is  the  JOTS  II  (Joint  Operational  Tactical  System)  barrier  model 
described  in  (4).  This  model  generalizes  slightly  the  computational 
method  in  Koopman  (2)  to  estimate  PCPfl(x),  the  probability  of  coming 
within  range  R  of  a  transitor  when  the  barrier  penetration  point  is  x. 

This  probability  is  given  as 


inhere 


PcprM  =  min  {fix),  1} 


(4) 


fix) 


2RK/T  for  |x|  «  T/2  -  RK 
(RK/T)  +  (T/2  -  |x |)/T)  for  T/2  -  RK  4  \x\  4 
(1/T)  V(R2~  (|x|  -  T/2)2)  [K2-])  for  T/2  +  R/K 
0  for  lx  I  T/2  +  R 


T/2  +R/K 
s<  |  x  |  N<  T/2  +  R 


K  =  Vi  +  (Vq/VtF  , 

and  x  is  0  at  the  center  of  the  searcher  track  T. 

If  the  barrier  penetration  point  is  given  by  a  probability  distribution, 
gx(x),  then 


(T/2)+R 

Pcpr  =  |  PcPfll^)  gx(x)  dx.  (41) 

-(T/2)  -R 

Currently,  the  SUBSUB  model  implicitly  assumes  that  the  barrier 
penetration  point  is  a  uniformly  distributed  random  variable  inith  mean 
0.  Using  (41)  allows  any  distribution  of  penetration  points.  For  example, 
an  attractive  strategy  for  the  target  would  be  to  attempt  to  maximize 
|x|  (that  is,  penetrate  as  far  as  possible  from  the  center  of  the  barrier), 
since  fix)  is  nonincreasing  in  |x|. 

tllhen  the  assumption  of  a  uniform  distribution  of  penetration  points 
with  mean  0  is  acceptable,  then  it  is  possible  to  closely  approximate  (41) 
with  a  closed  form  expression.  There  are  three  cases: 
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'  min  {2RK/T,  1}  for  BUU  «  T  -  2RK 
1  -  {(1  /BUU  T)  [max{.5(T+BUJ)  -  RK,  0}]2}  for  T-2RK  «BUU<  T+2R 
((T+2R)/BUU)  (1-  [rnax{T-R(K-1),  0}]2/(T2+2RT))  for  BUJ  »  T+2R 


P 


CPfl  - 


Luhere  K  is  as  before. 

Barrier  Search  Probability  of  Closure  (PrL).  The  SUBSUB  documentation 
describes  two  different  methods  for  estimating  PCL,  which  is  the 
probability  of  successfully  closing  a  target  detected  during  barrier 
search.  The  first  method  appears  to  compare  the  time  required  for  the 
searcher  to  approach  to  within  weapons  range  using  a  bearing  rider 
track  (i.e.,  the  relative  bearing  of  the  target  is  always  000°)  to  the  time 
required  for  the  target  to  escape  out  the  bottom  of  the  search  area. 

The  angle  R0BMflX  is  calculated,  which  is  the  maximum  angle-on-the-bow 
at  detection  allowing  closure  to  weapon  range  before  the  target 
escapes.  Then  based  on  a  geometrical  construction  which  is  equivalent 
to  assuming  a  cosine  distribution  of  target  ROB  at  detection,  the 
probability  that  the  actual  ROB  will  be  less  than  R0BMflX  is  determined. 
This  is  taken  to  be  PCL.  Rlthough  appearing  in  the  documentation,  this 
method  was  not  implemented  in  the  version  of  SUBSUB  received  for 
review. 

In  the  SUBSUB  code  (subroutine  BCLOSE),  a  more  direct  method  is  used 
to  determine  R0BMflX.  But  then  R0BMF)X  is  converted  to  PCL  in  a 
nonstandard  manner. 

The  program  assumes  that  the  searcher  requires  a  specified  length  of 
time  (variable  TIM)  to  conduct  target  motion  analysis  (TMR). 
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Immediately  after  detection,  the  searcher  heads  for  where  the  target 
will  be  at  the  end  of  the  TMfl  period.  (The  documentation  differs  from 
the  code  by  stating  that  the  searcher  closes  down  the  initial  target 
bearing.)  Rt  the  end  of  the  TMfl  period,  the  searcher  takes  a  corrected 
intercept  course  for  target  closure.  If  the  searcher  can  close  to  within 
weapon  range  before  the  target  exits  the  barrier  area,  then  the  initial 
ROB  at  detection  is  determined  to  be  small  enough  to  allow  closure.  The 
program  iterates  through  various  initial  ROB'S  (first  in  10°  then  in  1° 
increments)  until  the  maximum  ROB  allowing  closure  is  found. 

This  much  is  straightforward.  However,  the  program  then  converts 
R0BMrx  into  a  probability  of  closure  in  a  fashion  that  requires  further 
justification. 

One  "standard"  assumption  (which  is  implicitly  used  by  the  first 
procedure  described  in  the  documentation)  is  that  initial  ROB  will  have  a 
cosine  distribution.  Specifically,  reference  (3)  gives  the  probability 
density  function  of  ROB  at  detection  as 


f(0)  = 


.5  cos(9  -  y),  -tt/2  +  y  <  8  $  n  /2  +  y 
0  ,  otherwise 


where  y  -  tan_1(VC|/VT)  is  the  mean  ROB  at  initial  detection.  Given  a 
density  function  for  ROB  and  a  constant  flOBMF,x,  PCL  can  be  determined 
by  integrating  the  density  function  from  -R0BMF,x  to  RGBMF)X  (which  is 
particularly  simple  for  the  distribution  f(©)  given  above). 

In  SUBSUB,  PCL  is  computed  using  f(©)  when  the  searcher  detection 
range  is  sufficiently  small  (specifically,  less  than  T/(20  sin(R0BMflx)).  Arid 
for  larger  detection  ranges,  PCL  increases  linearly  to  a  maximum  value. 
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beyond  which  PCL  is  constant.  The  reasoning  behind  this  is  (apparently) 
that  the  assumption  of  a  cosine  density  for  initial  ROB  is  a  poor  one  for  ' 
large  enough  detection  ranges.  This  is  reasonable,  but  just  how  large  is 
"large  enough"  and  how  PCL  will  increase  requires  more  justification  than 
currently  appears  in  the  documentation. 

To  conclude,  the  subroutine  BCLOSE  in  SUBSUB  appears  to  give 
believable  values  for  PCL,  but  the  underlying  model  is  somewhat  ad  hoc. 
Further  analytical  justification  or  validation  with  simulation  results 
seems  appropriate  here. 

Prop  for  Rrea  Search.  SUBSUB  calculates  PCPfi  as  follows- 

Pcpfl  =  1  -  (1  -  (  jtR2/R))  exp(-2RVrt/(R-JTR2))  (6) 

where 

R  =  search  area  size, 

Vr  =  (V02  +  Vt2),/2  =  an  "effective"  relative  speed, 

R  =  detection  range. 

Problems  with  (6)  include: 

1.  Target  and  searcher  are  assumed  to  be  searching  in  the  same  area 
of  ocean.  This  seems  a  very  unreasonable  assumption. 

2.  Vr,  the  relative  speed  between  the  searcher  and  target,  is 
calculated  assuming  that  the  target  and  searcher  velocity  vectors  are 
perpendicular.  The  more  standard  assumption  is  that  the  angle  between 
the  velocity  vectors  is  uniformly  distributed  between  0  and  2rf.  Then  by 
the  law  of  cosines. 


2 If 

Vr  =  11/21T]  J(V02  +  VT2  -2VoVTcos0),/2de 

0 

tf 

=  (1/tt)  J  (V02  +  VT2  -  2V0Vt  COS0),/2  de  (7) 

o 

Evaluation  of  (7)  using  numerical  integration  might  not  be  advisable  on  a 
microcomputer.  However,  equation  (7)  is  can  be  approximated  by  the 
following  closed  form  expression: 

Vr  *  .38  rnax(V0,VT)  +  .62  (V02  +  VT2)1/2 

3.  Implicit  in  (6)  is  the  assumption  that  the  area  searched  the  instant 
the  search  begins  never  need  be  searched  again.  In  effect,  the  search 
area  reduces  in  size  from  fl  to  (fi-iTR2)  at  time  0+.  This  is  consistent  with 
a  target  stationary  in  relative  space  but  is  inconsistent  with  the 
assumption  of  a  moving  target  which  can  migrate  into  areas  previously 
searched. 

fin  alternative  expression  addressing  the  above  concerns  is 

Pcpfl(t)  =  11  -  (7fR2RST)/(BsRT))  exp(-2RVrRSTt/(fisflT)),  (8) 

where  Rs  is  the  size  of  the  searcher's  patrol  area,  RT  is  the  size  of  the 
target's  patrol  area,  and  flST  is  the  size  of  area  common  to  both  the 
searcher  and  the  target.  Rlso,  Vr  is  given  by  (7). 

It  is  noted  that  neither  (6)  nor  (8)  model  convergence  zones  Rs  a 
minimum  modification  for  convergence  zone  propagation,  the  irR2  term  in 
(8)  should  be  replaced  with  flC2,  which  is  the  area  of  ocean  giving 
convergence  zone  and  direct  path  detections.  For  example,  if 
detections  are  possible  at  ranges  of  0-4,  30-33  and  61-63  rirn  from  the 
searcher,  then 


1  K 


flcz  =  jt42  +  rf(332-302)  +  tT(632-612)  nm2 
=  1423  nm2 

Probability  density  function  for  tracking  range.  Equation  (1.25)  of  the 
SUBSUB  documentation  gives  the  density  function  for  tracking  range 

when  the  searcher  has  two  independent  sensors  as 

R 

<t>0(R)  =  <MR) +  -  ♦uitn)  ,f  ^R')  dR',  (9) 

-00 

where  <J>D1(R)  is  the  density  function  for  tracking  range  for  sensor  1  and 
4>D2(R)  is  the  density  for  the  tracking  range  of  sensor  2.  (Actually,  the 
documentation  gives  the  lower  Limit  of  integration  as  positive  infinity, 
but  this  is  assumed  to  be  a  typographical  error.)  Equation  (9)  is  not  a 
proper  density  function  since  it  integrates  to  a  value  greater  than  1 
(with  either  sign  of  the  lower  limit  of  integration).  The  proper  form  for 
(9)  is 

R  R 

4>o(R)  =  4>ot(R)  +  <t>D2(R)  -  4>m(R)  J^R')  dR'  -  <1>D2(R)  J>D11R')  dR'  (10) 

-00  -00 

Probability  of  Achieving  Track  given  no  Counterdetection  (TCP).  In  the 
SUBSUB  documentation,  TCD  is  calculated  by  conditioning  on  the  range 
required  to  achieve  "track"  on  the  target.  In  (1.27)  of  the 
documentation,  TCD  is  expressed  as 

J  (1-PcdIRD  <h(R)  dR 
o 

=  Prob{RCD  <  R track}-  HD 
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where  4>T(R)  is  the  density  function  for  tracking  range,  RCD  is 
counterdetection  range  and  RTrrCK is  tracking  range.  Given  the  event 
tree  structure  of  the  overall  program,  it  appears  as  if  the  correct 
probability  should  be  conditioned  on  the  searcher  achieving  a  secure 
detection.  That  is,TCD  should  be 

Prob{RCD  <  RtrrckI  rcd  <  Ro}-  02) 

Rs  long  as  the  target  has  some  secure  detection  capability  (i.e., 

Prob{RCD  <  Rd}  <  1),  equation  (12)  gives  a  value  strictly  greater  than  that 
given  by  equation  (11). 

fi  suggested  simplification  of  the  SUBSUB  model  would  be  to  assume 
that  TCD  is  some  user  specified  probability,  say  .8  or  .9,  depending  on 
the  searcher's  capability  to  conduct  passive  target  motion  analysis.  This 
seems  a  reasonable  assumption  and  would  reduce  the  number  of 
numerical  integrations  performed. 

Rttack  Probabilities.  Equation  (1.28)  of  the  SUBSUB  documentation  is  an 
expression  for  the  probability  of  a  secure  attack  (FCD).  This  probability 
is  calculated  by  conditioning  on  the  attacker's  launch  range.  The 
equation  is  correct  except  that  the  integration  should  be  performed 
from  minimum  launch  range  to  maximum  launch  range.  This  same 
apparent  error  occurs  in  Figure  1-6  of  the  documentation. 

Convergence  Zone  (CZ)  Environments.  CZ  environments  are  modelled  in 
a  nonstandard  fashion  in  SUBSUB.  It  does  appear,  however,  that  the 
final  effect  of  CZ  propagation  on  the  SUBSUB  calculations  might  be  small. 
In  CZ  environments,  PCPfl  is  modified  as  follows: 


PCpfi(R)  if  R<Rj 


(13) 


PcPfl(R)  = 


(PCPfl(R)  +  PCPft(Rcz}}/2  if  R^Rj 

Here  Rj  is  the  maximum  direct  path  detection  range.,  and  Rcz  is  the  range 
to  the  first  CZ.  The  justification  given  in  the  documentation  is  that  initial 
CZ  detections  will  make  subsequent  redetection  easier,  so  PCPfl  should 
be  larger  betuueen  direct  path  and  CZ  ranges.  The  result  of  (13)  is  to 
insert  a  discontinuity  in  PCPfl  at  range  Rr  See  Figure  5. 


Figures.  P  in  a  Convergence  Zone  Environment. 


This  modification  of  PCPfl  appears  odd  considering  that  in  non-CZ 
environments  PCPft  is  a  function  of  the  searcher  and  target  tracks  under 
the  assumption  that  detection  has  not  yet  occurred.  That  is,  PCPfi 
reflects  only  search  tactics  and  not  approach  tactics.  Rlso  the 
discontinuity  at  range  Rj  is  difficult  to  justify  intuitively.  Furthermore, 
(13)  apparently  applies  in  environments  with  only  one  CZ.  The 
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generalization  to  the  multiple-CZ  case  is  not  given  in  the  SUBSUB 
documentation. 

It  does  appear  that  the  effect  of  this  modification  of  PCFfl  might  have 
a  relatively  small  impact  on  SUBSUB  calculations,  at  least  in  the 
calculation  of  the  probability  of  a  secure  detection  (BCD).  This 
probability  is  given  by 

00 

Prob{RCDv<r}  fR[)(r)  PCpflU,r)  dr,  (14) 

o 

inhere  fR  (r)  is  the  density  function  for  detection  range  In  a  CZ 
environment,  fR[)(r)  mill  be  very  near  0  for  ranges  between  Rj  and  Rcz, 

since  most  initial  detections  will  occur  either  in  the  CZ  or  at  the  maximum 
direct  path  range.  But  it  is  for  those  ranges  between  Rj  and  RC2  that  PCPFi 
is  increased.  So  it  appears  that  the  CZ  modification  to  PCPfl  may,  in  fact, 
change  the  evaluation  of  (14)  only  slightly.  This  conclusion  might 
change,  however,  depending  on  how  the  multiple-CZ  modification  is 
accomplished. 

Analytical  Models  and  Discrete  Time  Step  Simulations 

SUBSUB  is  an  integrated  collection  of  analytical  models.  For  such 
models  to  determine  whether  or  not  a  particular  event  occurs,  an 
expression  for  the  probability  of  that  event  is  evaluated.  Then  a 
random  number,  uniformly  distributed  between  0  and  1,  is  drawn.  If  the 
random  number  is  less  than  or  equal  to  the  calculated  probability,  then 
the  event  is  said  to  occur.  This  is  opposed  to  discrete  time  step 
simulations  (e  g.,  SIM  II,  IBGTT  and  ENLUGS),  where  the  coordinates  of  the 
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individual  platforms  are  moved  during  each  time  step  and  the  status  of 
each  platform  (e.g.,  course,  speed,  depth,  radiated  and  self  noise  levels] 
is  updated,  fin  advantage  of  analytical  models  is  that,  in  addition  to 
determining  whether  or  not  the  event  in  question  occurs  this  time,  the 
calculated  probability  tells  the  user  houu  likely  an  occurrence  of  the 
event  will  be  next  time  and  all  subsequent  times.  This  could  be 
important  if  the  user  needs  to  know  whether  a  particular  outcome  is  a 
low  probability  event  or  something  to  be  expected  regularly. 

Rn  important  limitation  of  analytical  models  is  the  possible  lack  of 
verified  and  validated  submodels  to  calculate  the  required  probabilities. 
For  example,  in  SUBSUB  one  of  the  RSUJ  missions  modelled  is  systematic 
area  search  for  a  randomly  patrolling  target.  There  exists  no  simple, 
generally  accepted  analytical  model  for  this  search  scenario.  LUhat  is 
used  in  SUBSUB  is  a  modified  version  of  Koopman’s  random  search 
expression,  even  though  the  search  tracks  modelled  may  be  anything 
but  random.  It  this  a  reasonable  thing  to  do?  The  answer  depends  on 
the  particulars  of  the  situation  being  modelled.  If  the  relative  motion  is 
dominated  by  the  target,  then  random  search  is  probably  a  good 
assumption.  If,  on  the  other  hand,  the  target  speed  is  slow  relative  to 
the  searcher,  then  a  systematic  search  model  would  be  more 
appropriate. 

Or  consider  barrier  search.  The  model  used  in  SUBSUB  to  compute  the 
probability  of  detecting  a  transitor  assumes  (1)  a  constant  speed  search, 
(2)  the  target  track  is  orthogonal  to  the  searcher  track,  and  (3)  the 
barrier  penetration  point  is  uniformly  distributed  around  a  mean  position 
at  the  center  of  the  searcher  track.  If  it  is  desired  to  model  any  other 
case  then  SUBSUB  may  not  be  appropriate. 
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Or  if  environmental  conditions  are  changing  with  position  or  time  (such 
as  in  the  marginal  ice  zone),  then  SUBSUB  would  probably  not  be  the  best 
choice. 

Discrete  time  step  simulations,  on  the  other  hand,  are  more  robust 
and  can  generally  consider  more  types  of  scenarios.  If,  for  example,  it 
is  desired  to  know  if  a  transitor  successfully  penetrates  a  barrier  when 
the  searcher  conducts  a  sprint  and  drift  search,  the  physical  situation  is 
established  (initial  positions,  courses,  speed,  ranges,  and  depths)  and 
one  replication  of  the  simulation  is  run.  Either  detection  occurs  or  it 
does  not.  fln  advantage  of  this  method  is  that  the  modelling  can  often 
be  simple,  geometric,  and  straightforward  (compared  to  analytical 
models  which  frequently  become  rather  esoteric).  The  disadvantage  of 
this  method  is  that  it  does  not  give  a  probability  of  event  occurrence 
unless  many  replications  are  performed.  But  for  seminar  wargames, 
probabilities  are  often  of  secondary  importance.  LUhat  is  of  primary 
concern  is  typically  whether  or  not  a  particular  event  happens. 

To  conclude  this  section,  it  is  suggested  that  there  are  probably  some 
missions  and  scenarios  modelled  in  SUBSUB  which  could  be  handled  in  a 
more  straightforward  manner  with  a  discrete  time  step  simulation. 

Conclusions 

SUBSUB  is  a  large  program.  It  is  a  collection  of  analytical  models 
which  address  the  phases  of  one-on-one  submarine  engagements  from 
search  through  detection,  classification,  closure,  attack  and 
counterattack,  find  it  appears  to  run  reliably  on  an  IBM-PC 
microcomputer. 
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The  basic  structure  of  the  program  is  a  conditional  event  tree,  which 
is  hard  to  fault.  However,  some  of  the  component  models  do  have 
limitations  or  errors  which  should  be  addressed  in  the  documentation 
and  subsequent  revisions  of  the  program.  This  report  mentions  some  of 
these  problems,  but  there  may  be  others  which  escaped  notice. 

In  spite  of  the  program's  shortcomings,  it  is  judged  to  be  better  than 
the  microcomputer  programs  previously  used  at  NUJC  to  evaluate 
submarine  engagements.  It  is  recommended  that  the  current  version  of 
SUBSUB  be  used  by  NUJC,  but  that  the  known  errors  be  corrected  and  the 
technical  documentation  formalized  as  soon  as  feasible. 

One  separate  but  related  issue  briefly  raised  was  whether  analytical 
models  or  discrete  time  step  simulations  are  more  appropriate  for 
seminar  wargames.  There  is  probably  a  place  for  both.  If  analytical 
models  are  used,  however,  it  is  very  important  to  know  the  modelling 
assumptions  coded  into  the  expression  for  probability  of  event 
occurrence.  These  assumptions  might  be  more  restrictive  than  one 
would  like.  Compared  to  analytical  models,  simulations  tend  to  be  more 
robust  and  transparent  to  the  user.  However,  they  can  take  longer  to 
run  and  do  not  give  the  probability  of  event  occurrence. 
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